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Abstract 

We consider large fluctuations, namely overload of servers, in a 
network with dynamic routing of messages. The servers form a circle. 
The number of input flows is equal to the number of servers, the 
messages of any flow are distributed between two neighboring servers, 
upon its arrival a message is directed to the least loaded of these 
servers. Under the condition that at least two servers are overloaded 
the number of overloaded servers in such network depends on the rate 
of input flows. In particular there exists critical level of input rate that 
in case of higher rate most probable that all servers are overloaded. 

Keywords: Large Deviation Principle, Queuing Networks, Dynamic Routing. 



1 Introduction 

This work presents an effect in a network with interacting servers that can be 
called a spontaneous collective behavior in case of large fluctuations. 

We consider networks with dynamic routing of messages. In such networks the 
server to which a message is directed depends on the network's state at the message 
arrival moment. One of the problems arising here is the analysis of probability of 
large fluctuations, for example probability of large delays. 



There are many works where large fluctuations in networks with dynamic rout- 
ing have been investigated. In [I] - [9] the networks with two servers and three 
independent input flows have been considered where only one flow is divided be- 
tween two servers depending either on the workload of servers or on the queue 
lengths. In [TO] a network with a group of servers and several flows has been con- 
sidered where each flow is assigned to some subgroup of servers; upon its arrival 
a message selects a server with the shortest queue (i.e. a queue with least number 
of messages). In this work the large deviation principal for the flows upon the 
servers is proved. We would like to stress that the service time of messages in [10] 
is exponentially distributed. That allows to use Markov property for the flows to 
servers (after splitting of the input flows.) 

The more full list of references can be found in the mentioned works. 

Here we consider circle networks that are formed by k servers and k identical 
independent input Poisson flows. Messages of any flow are assigned to two nearby 
servers (see Figure 1) A message direction depends on the workloads on these 
two servers, namely upon its arrival a message is directed to the server with the 
smallest workload. The constant speed of work of each server is equal to 1, the 
discipline is "first in - first out" (FIFO). If a message finds the server busy it 
is put into a infinite buffer to wait for service. We consider the networks that 
work stationary. That means that with probability 1 the queues do not increase 
infinitely. But there may appear large fluctuations, for example during a short 
period one of flows may bring very large amount of work. We say that during this 
period the flow is overheated. Suppose for example that the flow f± is overheated. 
Then the buffers of servers s\ and Sk that are assigned to this flow will contain a 
large amount of work. What is the behavior of other flows ? 

In this paper we show that in case where f\ is overheated there exist at least 
two scenarios of network performance. What scenario is more probable and is 
realized depends on the rate of input flows. Namely, we show that there exists 
such value of input flow rate that in case when the arrival rates are above it the 
overheating of flow j\ coincides with overheating of all flows. Such behavior may 
be regarded as the spontaneous appearance of collective behaviour. On the other 
hand in case of low rate the overheating of f\ does not bring the overheating of 
other flows. There may exist also intermediate stages where the overheating of 
one flow coincides with the overheating of a number of neighboring flows. The 
existence of such intermediate stages depends on message length distribution of 
input flows. For some length distributions, for example for exponential one, the 
intermediate stages do not exist. 

The proof of this result is based on large deviation principle of [IT] (see also 
|12j). We do not present the details of this reduction. The details of application 
of large deviation principle for several other problems one can fined in (6] 
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Figure 1 



3 



|7]. Instead of the proof we present some ideas needed for the proof and some 
fragments of rigorous arguments. Therefore the following text is not a proof in the 
conventional sense. 

In the next section we introduce several concepts and definitions and describe 
the main ideas of our approach. In section 3 several examples are presented. 



2 Cyclic networks. Main result 

The cyclic network of size k, k > 3, is formed by k servers and k input flows. 
Let F = {fx, fk} be the set of flows and S = {s±, Sk} the set of servers. All 
servers work with speed 1, the discipline is FIFO. Each server has an infinite buffer. 
The messages of flow /j € F are served by two servers Sj € S and Sj_i G S (here 
always % =F 1 mod k). Each message of fi is directed to that of servers Sj, Sj_i 
which becomes idle first. 

Let the random sequence (£n) T n)) n = 1,0,1,..., describes the input flow 
fi £ F. The random variables are the intervals between the arrivals of mes- 
sages of flow f{. The random variables £^ are the lengths of messages. We consider 
homogeneous Poisson flows therefore all {Qii T n) are independent and equally dis- 
tributed, T % n are exponentially distributed with rate A, i.e. Pr(r^ > t) = e~ xt . The 
distributions of input flow do not depend on i. 

Introduce a random pair (£, r) that is distributed as any pair (£^,t^). Each 
sequence (C n > T n) 1S numbered so that Tq < and r\ > 0. 

We propose that £ has exponential moments, i.e. there exist such 6 + , < 
9 + < oo that 

<p{6) = Ee e€ < oo for 9 < 6+, 

lim (p(0) = oo. (I) 
0T0+ 

The condition that guarantees the stationarity and even ergodicity is 

\<f/(0) < 1. (2) 

This condition is intuitively obvious but its sufficiency for the existence of station- 
ary stage needs a proof. In particular one may construct a Liapunov function that 
shows that with large probability the workloads of all stations are located in a 
compact region. We omit this construction. 

To formulate the problem that is investigated we introduce a notion of virtual 
message. That is a message that arrives with flow fi at a nonrandom time moment, 
for example at and has zero length. It joins a queue at Sj or Sj_i, depending on 
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where the workload is less, and waits there its service. Let be the ith virtual 
message's waiting time. We are interested in the probability of large waiting for 
virtual massage that arrived with f±, i.e. Pr(oJi > d), where d is large. In some 
cases, for special choice of £ distribution the needed probability can be calculated. 
But in general case it is impossible to present the explicit expression. Therefore 
one looks for the asymptotics of probability 

Pr(wi > nd), (3) 

where n — > oo. More detailed , one looks for 

J = lim — lnPr(wi > nd). (4) 

n^oo n 

This problem belongs to the theory of large deviation. We use this theory even 
where it is possible to get an explicit solution. 



2.1 Collective fluctuations 

Next theorem states the existence of input flow values A& < X k such that separate 
different types of system performance. Let 



As we mentioned before the system performance is stationary if A < A. 

We say that the flow /j is overheated at t = if Tj(n) = (u>i > nd). It means 
that a virtual message that arrived with fi is waiting for being served at least nd. 
Let Ti(n) be the complement to event Tj(n). 

Theorem 1 . For any system size k, k > 3 

and for any d there exist A&, A fe (not depending on d), such that < A& < X k < 
A, and 

► if X < Xk then 

1) Pr (Ti(n) n* =2 Ti(n) \Ti(n)j -»■ 1, as n ->■ oo. 

2; J = lim n _ +00 ^lnPr(r 1 (n)n- =2 r i (n)) = 20(A, l)d 
where 9(X, 1) is a positive root to equation 

20 = \[<p(9)-l], (5) 



5 



► if X> X k then 

3) Pr (nf=i Ti(n) ^(n)) -^losn^co, 

4) J = Hm n _ >00 =± In Pr (f|- =1 ^(n)) = fcfl* 
where 6* is a positive root to equation 

9 = \[<p(0)-l]. (6) 

It follows form this statement that in case A < A^ the main contribution to the 
probability of event Ti(n) = (u>i > nd) brings the flow fi which the virtual message 
arrives with. Only this flow is overheated, the others stay not overheated. But 
in case A > X k all input flows are overheated. Though the virtual message is 
combined with flow f\ only its large delay is connected with an effect of collective 
behavior that is similar to the effect of spontaneous magnetization in statistical 
physics systems. 

It is possible to present a more accurate statement. In the next section we 
introduce several definitions and present a precise theorem. 

2.2 Random processes connected with the system and 
large deviations 

Below we introduce several random processes connected with the cyclic network. 
1. Process 

at) = E & 

describes the amount of work brought by fi during time interval [0, i). 
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2. k dimensional Markov process w 1 , ...,w k denned by generator 

Lg(w (1 \...,w^) = 

= -E5)V>0)(» (1) ) + 
i=i 

k 

+ XY,[mw\...,w' l + ^...,w k )- 
i=i 

-g(w\ w\ ...,w k ) I^i^i+i^w 1 , w i+l ) + 

k 

+ \^2[Eg(w^,...,w l+1 + t-,w k )- 
i=i 

-g(w\ w i+1 , ...,w k ) I (tt , i>u;i+ i)K, w m ), 

describes the loading of all servers. We propose that w l {0) = 0, i = 1, fc. 
The symbol is the indicator of a set A, § is arbitrary differentiate function 
in L2(M fe ). The second and third sums of the right hand side indicate that 
a message of length £ arriving with is directed to one of servers Sj, Sj+i 
depending on where the loads w ! , tt; l+1 is less. 

A typical trajectory of this process is a k dimensional function with non 
negative components. Each component is formed by jumps and piecewise 
linear functions between the jumps. On the intervals where the trajectory 
is differentiable its derivative is either — 1 (if it is positive) or (if it equals 
0). 

3. Process (w 1 ^), ...,w k (t)) describes the amount of work that was brought to 
each of k servers presuming that w l (0) = 0. Let t l r > 0, r = 1, ... be the set 
of jumps of process w % and a*, r = 1, ... the sizes of these jumps. Then 

w\t) = £ 4. (7) 

r: 4<t 

The trajectories of this process are step functions, where jumps coincide in 
size and time with those of process w 1 , w k . 

We define these three fc-dimensional processes ( = £ fc ), w = (w 1 , ...,w k ) 

and w = (w 1 , ...,w k ) onto the same probability space. Here the jumps of w and 
w coincide in size and time. Let Q % have a jump at time moment t. Then at the 
same moment either u> l_1 or w l has a jump of the same size, depending on whether 
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an inequality — 0) < w l (t — 0) or the opposite one takes place. In case of 

equality we propose that process has a jump. 

Defining the processes on one probability space one can define a mapping 

G : X° -> X° (8) 

of the set realization of processes £ onto the set realization of processes w. Here 
X° is the set of nondecreasing stepwise functions that are equal to at 0. 

All event that are considered here are connected with asymptotic characteris- 
tics of network as n —* oo. Therefore we introduce the scaled versions of processes 
and events. 

We start with introduction of scaled versions of the processes [T] and [3} 
1. 

C(t) = -C(nt). 
n 

2. 

= l w\nt). 
n 

Being interested in asymptotic behavior of probabilities of different events we 
introduce also several notions defined by the sequences of events. 
Let {xT, a (-)} be a piecewise linear trajectories 

*r.(i) = K iftG[ °' T] 0) 
' |V(0) if t >T V 7 

We say that a flow /, is overheated on an interval [0, T] if following event takes 
place 



■Aln(x) = SUp \x(t) — (n(t)\ < £ ) , X = XTa, O>l,e>0. 
V* e [°' T ] / 

Remark that the notion of overheated flow is introduced only for piecewise linear 
trajectories. We do not need more general definition. 

We say that a flow /j is not overheated on an interval (0, T) in case of event 

"4e,nO) = ( su Pte[o,T] \ x {t) - C(*)| < e ) i x = x T,a, a < 1, e > 0. If a = a (0) = 
X(p' (0) then not overheated input flow is close to its mean value. 

Remark that the distribution of uj\ coincides with the distribution of 
min {sup t>0 { , u? 1 (t) — t}, sup t>0 {w k (t) — t}}. That can be shown easily, we omit 
the needed construction. (For the problems considered in |6l [7] similar facts 
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are explained in these works, the needed construction for a one-channel systems 
one can also find in monograph [15] ). 

The probability (|| can be expressed in terms of scaled processes w n (t). It is 
equal to 

Pr(wi > nd) = Pr (M n ) , (10) 

where event 

M n = ^min jsup ju^(t) - t| , sup -t}| > d 

means that both processes u>k intersect the line u(t) = d + t. 

Using the mapping G of process ( n realizations onto process w n realizations 
(see (18])) we can introdus an event M n = G~ 1 M n - 



Theorem 2 . For any size k, k > 3 of the network and for any d > there exist 
not depending on d Xk and X k , < Xk < X k < X, such that 



► if X < Xk then there exist and such that 

1) lim e ^o lim n _ Pr (^(^.a; ) fliU -^.nO^a 



k 



1 



2) J = lim £ ^o hm n ^o In Pr [A l ^ n {x T ^ a * ) f| i=2 ^l, n ( x T*,a(o) ] 
20(A,l)d, 

where 9(X, 1) is the positive root to equation 



3) besides 



20 = Afo>(0)-l], 



Ay/(0(A,l))-l' 



► if X > X k then there exist T£ and a* k such that 

M n 

5) J = lim e ^o lim n ^ ^ In Pr (f) k =1 Al >n {x T *, a 
where 9* is the positive root to equation 



4) lim e ^ hnin^o Pr ( fli=i A\ 



k9*d, 



(11) 
(12) 



X[cp(9)-1], 



(13) 
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6) besides 



(14) 



This more precise form of the theorem presents the mean dynamic of con- 
ditional process under the condition of event M n - Namely, the trajectory 
(xT*, a * , a(°H, ...,a,(°H) is the conditional mean dynamic of processes £*) on 

the interval [0, T*] in case A < Afc, and the trajectory (xt* ai>— > X T?, at) is the 
conditional mean dynamic on the interval [0, T£\ in case A > \ k 

2.3 Ideas of proof 

The value of J could be found if the large deviation principle for processes u; 1 (i) 
and w k (t) would be known. These processes are not Poisson. But they are the 
functionals of C l {t), ...C, k {t). We remind that they are defined on the same proba- 
bility space. The processes i = 1, k are Poisson. Therefore one can use some 
known results on large deviations of Poisson processes. We use the results on large 
deviation principle from [11] (see also |12]). 

For the problems that are considered here the application of large deviation 
principle consists of two parts. First one has to check the validity of the large 
deviation principle and describe the corresponding to the problem event. After 
that one has to find minimum of rate function on the event. As usual, application of 
large deviation principle reduces the problem to the search of the point (sometimes 
several points) of the event where the rate function take minimal values. A small 
neighborhood of the point of the minimum brings the main contribution to the 
asymptotic of logarithm of event probability. In our case the event is some set 
of trajectories and the rate function is an integral functional on the trajectories. 
Therefore a search of trajectory that minimizes the rate function is reduced to a 
variational problem. 

To check the validity of large deviation principle one has to present the topo- 
logical space, the sequence of measures on this space that tend to a 5-measure and 
the rate function. 

The topological space is a set of non decreasing k dimensional functions on 

[o,«o 

x = x k = {(xH-),..,x k (-))}, 

equipped with uniformly week topology. 

The detailed description of these functions one can find in |llj . where the 
uniformly week topology is introduced on this space. We use this topology. It is 
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explained in the same paper why the week topology is too week and one has to 
use the uniformly week topology for the problems of the considered kind. 

The sequence of measures P n is defined by processes ( n . 

Finally, the rate function / is defined on (X, P n ) by equalities 

k 

I{x\...,x k ) = Y,I^), (15) 

i=l 

where 

/•oo 

I(x l ) = / sup {9x\t) - \[<p(6) - l]}dt, (16) 
Jo e<e+ 

if x l are the absolutely continuous functions. We do not give a definition of func- 
tional on not absolutely continuous functions because condition ([TJ permits to 
avoid them. 

The event T where minimum of rate function we have to find is defined by the 
variational problem presented below. 

First let us extend the mapping ([8| onto the whole set X of not decreasing 
functions 

G : X -► X. (17) 

The extension is denoted by the same symbol. 

Consider X = (x 1 ,...,x k ) S X and the functions a l (t) G [0,1], i = l,...,k. Let 

y \t) = a l (t)x\t) - (1 - a l+1 (t))x i+1 (t), i = 1, k. 

Here Y = GX = (y 1 , ■■■,y k ) is a solution to the following optimization problem 

E = mm^2 J \y\t)-y^\t)\dt, (18) 

i 

a = (a^t),...,^*)). 
The event is 

T = \X = (x 1 ,...,x k ) : mm\sup{y\t) -t}, sup {y k (t) - t] 1 > d, Y = GX) . 
I U>0 t>0 J J 

(19) 

It follows from the routing rules that the flows {y l ) upon the servers given the 



input flows (x l ) are defined by solution to (18). That is the consequence of the 
system routing rules. We omit the proof of this fact. 

To prove the theorem one needs to find the minimum 

= min {/(X) : Iff}, (20) 
11 



As the rate function is an integral functional on X, the solution to (20) is re- 
duced to the solution of variational problem on the set of non-decreasing functions 
from T . 

In fact it is sufficient to find ( |20[) on the piecewise linear functions of form 
Q. That is because: 1) condition ([I]) permits to restrict oneself by absolutely 
continuous functions; and 2) for the homogeneous processes with independent 
increments the rate function on the trajectories that connect two fixed points 
during a given period takes its minimum on the line that connects these points. 
We omit the detailed explanation of these facts (see [HI [13]). 

Below we use the following not formal terminology. The trajectory X = 
(x 1 ,...,x k ) is called an input flow and the image Y = GX = (y 1 ,...,y k ) a load 
flow . A small neighborhood of trajectory X contents "the real" (jump-wise) tra- 
jectories of input processes and a small neighborhood of trajectory Y contents 
"the real" trajectories of load processes w l n . 

A trajectory X = (x 1 , ...,x k ) G X that represents k input flows will be called 
a input configuration or simpler a configuration We consider only trajectories of 
form ([9]), each is characterized by a pair (a l ,T l ) where a 1 = x l (t) on the interval 
[0,T*]. We denote by X = {a, T} a configuration where all T % are equal, here 



a = (a 1 , a k ) 



For functions of form ([9| a solution to ( 18 ) is equivalent to solution to following 
optimization problem 



Ol D = min^|6 i 

where 



a 

i=l 



i) i = aV + (l-a w )o i+1 . 
a* G [0,1]. 

a = (a 1 , a k ) {i + 1 mod k). 
Obviously the solution exists. 

In our problem for given input flow configuration X = {a, T} vector b repre- 
sents the load configuration Y = {b, T} . 

In addition to the above correspondence of all X with Y we consider also the 
correspondence of subset of input flows with a subset of load flows to the assigned 
servers. 

We are interested only in connected subsets F = {/ r +i, fr+l} ^ F of in- 
put flows. Let X = {x r+1 , x r+l } be a configuration of such flows. It is sup- 
posed that x l G X are of form all T are equal, i.e. X = {a r+ i^ r+ i,T} 
where a r+ i^ r+ i = (a r+1 , a r+l ) G and a 1 = x l (t),t < T. The corresponding 
configuration of load flows to the assigned servers S(X ) = {s r , s r +i, s r+ i} is 
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Y' = {6 r>r+ ,,T} = {{b r ,b r+1 ,...,b r+l ),T}. Here vector b r , r+l G R l+1 is defined by 
a solution to optimization problem 

r+l 

D r , r +i = J™ 11 ^ — b 1 

l=T + l 

where 

'(l-a i+1 )a m if i = r, 

(1 - a i+1 )a i+1 + aV ifr<i<r + Z-l, 

aV if i = r + Z, 

e [0,1], 

a rjr+i = (a r , ...,a r+ '). 

We look now for solutions to 02 in case D r ^ r+ i = 0. Below the notation D x F > 

is used to indicate that D is calculated with respect to a set of input flows F' 
for a given configuration of all flows X. A connected set F is called balanced if 
=0. 

As we are interested in large deviations keeping in mind f\ only subsets F C F 
with fi G F are considered. 

A balanced F CFis said to be maximal if D x p n > for any F D F . A 
configuration X = {a, T} may posses several maximal balanced subsets. 

Our aim is to find for a configuration I = {a,T} £ f for which the rate 
function is minimal. 

The rate function for X = {a, T} is equal to a sum 

k 

I(X) = T^(0 i a i -\[p(e i )-l]), (21) 
i=l 

where 0j are defined by a 1 = \<p'(0i). 

Let X = {a, T} be a configuration, F a subset of input flows and Z a set of 
configurations such that 

Z(F') = {Z = {c, T}: c l = a\ if h G F'}. 

Rate function for this set is 

I(Z(F')) = I(X ) = T (^-A[^)-l])- (22) 

where Ao = {a, T} is a configuration with a l = a{ if /i G F' and a J = a< ) if 
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For balanced with respect to X = {a, T} subset F' we consider 



E 



i a 



°X,F 



\F'\ 



Here \F | is the number of flows in F . 

Let h > 1 and let Z(h,F ) = {Z = (a, T), h ZF > = h} be a set of trajectories 
with fixed mean value h on F' . 



Lemma 1 For a set of configurations Z(h,F) a rate function I(Z(h, F )) on 
interval [0, T] is equal to 



(23) 



I{Z{h, F)=(h0- \[<p(6) - 1] ) T\F 
where is determined by equality 

h = \(p'(6). 
If F' is a maximal balanced set then 



lim lim —lnPi(A £n {h,F) \ =I{Z{h,F')). 

£—►0 n— too n 



(24) 



(25) 



Here we denote 



A £ , n (h,F')= |J 

ZeZ(h,F') 

Z = {c,T} = {z\...,z k }. 



f^l -^£,71 ( z T,c l ) ^ f^l -^£,n( z ^ ja (0)) 

.^/^eF' r-fjeF\F' 



This Lemma shows that equal overheating in all flows from F' is "more" prob- 
able than the not equal one. We write "more" because that is an asymptotic 
result. In fact the considered probabilities decay exponentially and the exponent 
determined by rate function / is minimal in case of equal flows. 



Proof. Rewrite the expression for (21) 



I(Z(h,F'))= {(h + g t )e i -X[i P (0 l )-l])T, 

i-.fiEF' 

where g% = c l — h. The solution to system = is: Q{ = 9j for all i,j : 

fi, fi G F' . The solution is unique thanks to monotonicity in 8 of ^( c6> ~ A ('^( e )~ 1 )) ; 



c = \ip'(9). Equality (25) followes from (22). 
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▲ 



It follows from X £ J 7 that 



bT-T>d, 

where b = min{b k , b 1 } and b is a solution to Ol. 

Introduce now a set of configurations C(h,F') where F' is a single maximal 
balanced set for Z = {c, T} G C(h, F'), c { = h, G F' , h > 1. Let 



c= U UUc(",n 

F 1 <ZFh>lT>0 



(26) 



We want to find the rate function for C. To find 1(C) as A is fixed one has to 
minimize C in T, h and F (see (23), (24)). Because of circular structure of the 
network instead of minimization in F one can minimize in number I = \F \ of 
input flows F . If S' is the set of servers assigned to F' then \S'\ = I + 1 as I < k 
and \S'\ = k as / = k. 

Below the calculations are based on the following argument: the overheat of 
flows F' that bring the overload of servers S' can be considered in case I < k as 
a overload of a one-channel system with a server of speed I + 1 and input flow of 
rate XL In case I = k a one-channel system has a server speed k and the input 
flow rate Xk. We call such a one-channel system an auxiliary system. 

Consider C(F') = \J h>1 \J T>Q C(h, F') and look for I(C(F')) as \F'\ is fixed. 

Suppose first that I = \F \ < k — 1. Obviously F is assigned to \S \ = I + 1 
servers. We can consider only such load flows that form configuration Y = {b, T}, 
with b % = b = A^h, Si G S . This configuration and its small neighborhood belongs 
to the manifold T if 

-hT - T > d 



Now by (23) and (24) we get 



I(C(h,F 



inf 

T 



l + l 



9 ti 



\[<p{0) IT 



(27) 



where 9 is defined by ^ — ^ = A<^9 (#). It is easy to see that (27) has its infinum 

as 

(l + l)6 = lX[<p(6)-l}. (28) 
Denote by #(A, I) a positive solution to the last equation. Then we get that 

J(A, I) = inf I{C{F))) = {l + 1)0(A, i)d (29) 
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and the optimal b and T are 

I 



6 = 6(0 = A— <p'(9(X,l)), T = T(l) 



X^'(9(XJ))-1 



(30) 



In case where I = k the sum of servers speed is k and the sum of input flow rates 
is Xk. The configuration Y = {b, T} of flows to the servers is such that b l = b = h 



for S{ 6 S and we have the inequality hT — T > d. Therefore (27) becomes 

T+ --X[<p(9)-1] 



I(C(h,F)) =Tk[6 



Optimization in T and h gives 



J(X,k) =inf I{C(h,F)) = kd*(X)d, 



where 6*(X) is a positive root to 
The optimal b and T are 



\[<P0)-1\. 



b = b* = X(p'(9*(X)), T = T* 



d 



Xif'(9*(X)) - 1 



(31) 
(32) 
(33) 



Remark that J(A, k) < J(A, k — 1) for any A. Really, j£\L\ = I = k — 1 then all 
servers are overloaded, J(A, k — 1) = k9(X, k)d, where 9(X, k — 1) is a solution to 
9 = X^[ip(9) - 1]. At the same time J(A, k) = k9*(X)d, where 0*(X) is a solution 



to @. Therefore 6»*(A) < 6>(A, — 1), see Fig.2. 

To get rate function 1(C) as A; and A are fixed we have to find I, I < k, that 
brings 

min 



k9*(X), mm(l + l)9(X,l) 

l<l<k 



X(<p(9) - 1) 




A 2 > Ai 



Fig. 2. 
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Lemma 2 1) For any k, I, 1 < I < k — 1 there exist such Xf.,, < X kl < A 
i/iai J(A, /) < J(A, fc) as A < X* kl and J(A, /) > J(A, A;) as A > X* kl . 

2) For any l\, I2, 1 < h < I2 there exist such Xi 2 / 1; < A; 2 ; 1 < — A i/iai 
J(X,h) < J(X,l2) as X < A; 2j / 1 and J(A,Zi) > J{X,l2) as X > A/ 2 ^ 1 . 



Proof. 1). For the start we find such A = X* k b I < k-1 that J(X,l) = J(A, k). 
Denote (I + 1)9 = $1, I < k, k9 = It follows from ([32]) and (pi) that 



r k = xk 



n 

A: 



1 



^ = /A 



til 



l + l 



1 



(34) 



By (29) and (31) the needed equality is achieved as = Let us show that 



equation 



l + l 



1 



k 



has a unique solution ; and 

< K,i 

Really, 93(0) = 1, therefore 



i) 



KMJk) - 1; 



1 



< A. 



I 

lim - 
k 



l + l 



< 1. 



Further, if 9 + < 00 then 



lim (p 



i) 



$/(i+i)^e+ \l + l 
And if 9 + = 00 then 



00, lim (p. 

\k 



< 00. 



lim 

1?— >oo 



nf 



00. 



(35) 



That means that there exists $ = $%i such that (35) takes place. The uniqueness 
follows from because (p and all its derivatives are convex.. 

The function 6>(A), < 9 < 9 + presented by (32) is defined for < A < A and 
9(X) monotonically decreases in A, (see Fig. 2). Therefore there exists such A^, 
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A£ £ < A, that corresponds to i? = $? > and for which the conditions of Lemma are 
fulfilled. 

2) . The proof of this item follows once again from existence and uniqueness of 
solution to 



h + l 



1 



Zv 



/ 2 + l 



1 



therefore it repeats the proof of item 1). But here we have to notice that 9(X,l), 

l + l 



< 9 < 9 + , presented by (28) is defined for A < A— — (in (32) we had A < A). 
Thus < Ai 2) j a < X 1 - 2 ^ and in general case it may be that A^^ > A (see Fig. 3). 



The proof of Theorems 1 and 2 follows from the Lemmas. 
Let us set 

A fc = min[A^,A Zjl ], 

Kk 



X k = max A 
i 



By Lemma 2 J(A, 1) < min J(A, /) as A < min^ A^ i and therefore for fixed k we 

l<Z<fc 

have J(A, 1) < min 1< ; </ t{ J(A, I), J(A, k)} as A < A&. The way we get (29) indicates 
that as A < Afc the statements 1) and 2) of the Theorem take place and the values 
(30) correspond to the values of idem 3) of the Theorem. 

Further, J(X,k) < min J(A, /) as A > X k . The way we get (31) indicates that 

l<Z<fc 

as A > X k the statements 4) and 5) of the Theorem 2 take place and the values 
(33) coincide with the values of item 6) of the Theorem 2. ▲ 



3 Different distributions of message length. 
Examples 

It is clear that if Afc = X k for some k then as f\ is overheated then, depending on 
A, most probably either only fi or all flows are overheated. The questions are: 
when Afc = X k , what happens if A& < X k ? 



To answer these questions we look at the location of curves $i and i9f (see (34 )) 
on (A, #) plane. 

Preposition 1 If A/ 1 > A for any 1 < I < k, then Xk = X k , X k — > A as k —>■ oo. 
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Proof. By Lemma 2 each pair of curves f?i(A), i?z(A) defined by, (34) 
has a unique point of intersection (A^i, #^1) and i?i < Z > 1, as A < A/^. 
Therefore from A^i > A follows that J(A, 1) < J(A, Z) for all A < A. Thus 
min J(A,Z) = min[J(A, 1), J(A, fc)], i.e. Xk = A fc . The curves i?z(A) do not de- 

l<Z<fc 

pend on /c, the curves $jjl(A) increase with fc, therefore A£ x = A fc — ► A (see Fig. 3). 




Fig. 3 































A 


2A 




b) 





Preposition 2 If $i(\o) > $i(\o) for some Z, Z > 1 and Ao < A then Xk < X k 
for sufficiently large k and X k — > A as k — > oo. 

Proof. It follows from the condition of Lemma that A; i < A. Remember 
that i9£(A) increases in fc. If Zc is sufficiently large then the pairs #i and 
i?z intersect at points A^ and A^ where A£ z > A^i, A^ 1 > A^i, and X kl > 
X* kv Therefore i?*(A) > #i(A) > 0j(A) as A M < A < X* kl (see Fig. 4),' i.e. 
mini<z<fc J(A, Z) < min[J(A, 1), J(X, k] for such A and that means A^ < X k . As in 
Preposition 1 A fc ^ A as A; ^ oo. ▲ 
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1. Exponential distribution of message length. 



The density of message length distribution is ce cx , c > and (f(6) 



the condition of stability is A < A = c. The equations (34) are linear having the 
form 

■& l = c + l(c - X), #% = k(c-\). 

It is clear that $i(A) < i9/(A) as A < c = A. On (A, i?) plain all $i(X) intersect at a 
point (c, c) and X k = X k = c|^f (See Fig 3 ). 

2. The density of message length distribution is 

1 



{c + g)e 



-{c+g)x 



+ (c- g)e 



(c-g)x 



, g < c. 



0(c — 9) c 2 — g^ 

Here we have — 1 = -, ^ ^ • The stability condition is A < A = . 

(c — 6) z — g 

( 

u-ef 



For A < A we can find tfj(A) = (I + 1)6>(A) using (|28| 

' A(c-0)- 5 2 
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the needed expression is 



0i(A) = (Z + l) 



XI 



X 2 l 2 



2(1 + 1) 



4(1 + 1) 



+ 9 2 



c + l 



X 
2 



A 2 



+ 9 



2 ([+11! 



(36) 



The value in square brackets of (36) increases in I, thus $j(A) increases in / as 
A < A. That means that all i?z(A) intersect as A > A (see Fig 3 b). By Preposition 
1 that indicates that A^ = X k . 

This example demonstrates the possibility of the following scenario: the over- 
heat of I connected flows may bring not only overload of I + 1 assigned servers, but 
also overload of other servers by not overheated flows. 

For example, let k = 3, c = 1, g = 0.5, here A = 0.75. The numerical 
estimation gives A3 1 ~ 0.418, i.e. 2Ag 1 > X. That means that as A is inside an 

interval f |, A3 -A then the overheated /1 brings not only the overload of s±, S3 but 
also the overload of S2 fed by not overheated /2 and f% . And it is easy to estimate 
that si, S3 are overloaded with greater speed than S2- 

3. Constant message length. 

Let Pr(£ = x) = 5(1/ 'c), c > 0. Here <f(0) = e 9 ^ c . The stability condition is 
A < c. 

c we get by (34) that = l( e ^K l+ ^ - 1). It is easy to 



For A : 
check that 



= A = 
01 (c) 



2, i.e. i?i(A) > 0z(A) as I is sufficiently 



> 2, and lim 

C i^oo C 

large. By Preposition 2 A& < X k as k is sufficiently large. 

The numerical estimates performed for c = 1 and k < 35 indicate that the 
behavior of A& and X k changes as k increases (and so do also the scenario of 1?; and 
$* k intersection). 

a) Here A& = X k = A^ x as 3 < k < 12, Xk increases in k. 
For example A^ ~ 0.311 as k = 3; A^ ~ 0.667 as k = 5; A& ~ 
X k ~ 0.883 as k = 12. 



0.857 as k = 10; 



b) As k > 12 the value A& does not change and is equal to A2,i ~ 
increases in k, Xk < X h . 



0.888; X k 



c) Up to k = 28 min J(A, I) = J(X, 2) as A^ < A < A fe . Therefore most probab 

two flows are overheated. For example X k ~ 0.910 as k = 15; X k ~ 0.935 as 
k = 20; X k ~ 0.940 as /c = 25. 
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d) As k > 29 in addition to interval (Xk,^3,2), where min/ J(X,l) = J(A, 2), 
there appears an interval (A3 5 2, X k ) where min J(A, /) = J(A, 3), and most probable 

three flows are overheated. Here A 3j2 ~ 0.956; X k ~ 0.959 as k = 30; A fc ~ 0.965 



These estimates show that nonzero interval (Afc, A fc ) is small and X k is close to 
A as k > 12. Presented numerical data and some analytic investigation suggest 

that as k — ► oo a "jump" from I to k (where I = arg min J(m, A) changes its 



All numerical data are presented with accuracy 0.0005 

We want to remark that in case of large fluctuations the collective behavior of 
dependent servers may take place for others, not circular networks. 
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